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ABSTRACT 



Maps of the peculiar velocity field derived from distance relations are 
affected by Malmquist type bias and selection effects. Because of the 
large number of interdependent effects, they are in most cases difficult 
to treat analytically. Monte Carlo simulations are used to understand 
and evaluate these effects. In these simulations the "true" spatial dis- 
tribution and relevant properties of galaxies as well as selection effects 
and observational uncertainties are realistically modeled. The results 
of the simulation can directly be applied to correct observed peculiar 
velocity maps. The simulation is used to investigate biases in samples 
of measured peculiar velocities by Lynden-Bell et al. (1988), Willick 
(1990) and the new sample of spiral galaxies by Haynes et al. (1993). 
Based on the results obtained from the application of our method to 
toy models we find that the method is a useful tool to estimate the bias 
induced both by inhomogeneities and selection effects. This is a crucial 
step for the analysis of the Haynes et al. sample which was selected with 
a redshift dependent criterion. 

Subject headings: large scale structure of the universe - cosmology: ob- 
servations - galaxies: clustering - galaxies: distances and redshifts 
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1. Introduction 

Maps of the peculiar velocity field of galax- 
ies are a powerful tool to test scenarios for 
formation of large scale structure. However, 
the interpretation of such maps must take into 
account biases introduced by the coupling of 
errors in the distance measurements and the 
spatial distribution of galaxies and/or selec- 
tion effects in the samples considered. Such 
biases in the derived distances are present 
even if the underlying distance relation is ex- 
actly known (which we will assume through- 
out this paper). A discussion of how a bias- 
free distance relation can be derived from 
data will be presented by Giovanelli et al. 
(1995). Biases originating from the spatial 
distribution are usually referred to as Malmquist 
bias and to correct for it requires a detailed 
knowledge of the underlying distribution of 
galaxies from which the sample is drawn. Bi- 
ases introduced by the selection of the sample, 
which we will call selection biases through- 
out this paper, are caused by the correlation 
of the parameters used to select the sample 
with parameters used to estimate distances. 
In order to correct for such biases, the dis- 
tribution function describing relevant galaxy 
properties is needed, as well as the detailed 
criteria adopted in selecting the sample being 
considered. 

A number of correction schemes for Malmquist 
bias have been suggested by different workers 
{e.g., Feast 1987, Lynden-Bell et al. 1988; 
Landy & Szalay 1992; WiUick 1994). A com- 
mon feature of all these approaches is that 
they require a priori knowledge of the dis- 
tribution of galaxies in space, which has to 
be either assumed (e.g. Lynden-Bell et al. 
1988) or, in the non-trivial case of an in- 
homogeneous distribution, derived from the 
data before the correction for the bias can 



be made. Different methods for estimating 
the required density distribution have been 
considered. Dekel (1993) built a treatment of 
the inhomogeneous Malmquist bias into their 
reconstruction of the density field, assuming 
that the distribution of galaxies traces the 
underlying matter distribution derived from 
peculiar velocities. Another approach uti- 
lizes the self-consistently reconstructed den- 
sity field derived from redshift surveys {e.g.., 
Dekel 1994, Hudson 1994). 

Biases resulting from selection effects can 
in some situations be either comparable or 
larger than the Malmquist bias. To correct 
for selection biases, detailed knowledge of the 
selection process is necessary. However, even 
if this information is available it is difficult to 
derive simple expressions for situations where 
the parameter used to select the sample is not 
identical but correlates to one or several pa- 
rameters used to derive distances. This is es- 
pecially true when a complex set of different 
selection criteria is used to select the sample. 
For example, samples where the selection de- 
pends explicitly on the redshift, like the re- 
cently completed survey of peculiar velocities 
of Sc galaxies (Haynes et al. 1993, hereafter 
referred to as the the Sc sample), arc partic- 
ularly difficult to correct for this kind of bias. 

The biases mentioned above arc strongly 
related to the scatter of the particular dis- 
tance relation used. For example, the scatter 
in the Tully-Fisher (TF) relation depends on 
the line widths of the galaxies in the sense 
that the scatter is large for galaxies of lower 
luminosity (Giovanelli et al. 1995). Such an 
effect would be difficult to treat analytically. 
This fact and the unique selection criteria of 
the Sc sample have motivated us to address 
the question of biases, both Malmquist and 
selection biases, using a Monte-Carlo (MC) 
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technique. 

MC simulations of Malmquist bias and/or 
selection effects have been used in the past 
by a number of workers (e.g. Landy & Szalay 
1992, Teerikorpi 1993, Wilhck 1994), mainly 
for the purpose of testing analytic corrections. 
In contrast, here the MC technique is used 
to determine the bias at each point in space 
to correct measured peculiar velocity fields. 
The essence of our method is to simulate 
galaxy catalogs using the spatial and para- 
metric distribution derived from observations, 
and select redshift-distance samples similar 
to those currently available. The strength of 
the method is that all relevant effects can be 
incorporated to predict the bias without the 
need for any simplifying assumptions. 

In section 2 we review analytical correc- 
tions for Malmquist bias which are used to 
test our MC results in simple cases. In sec- 
tion 3 a detailed description of the technique 
is presented. In section 4 we compare the bias 
derived from our MC simulations with analyt- 
ical corrections and in section 5 we derive the 
bias field for observational samples. A brief 
summary is presented in section 6. 

2. Biases in the Peculiar Velocity Field 

Different analytical expressions to correct 
for the Malmquist and selection type bias can 
be found in the literature (see e.g. Willick 
1994; and references therein). These expres- 
sions depend critically on the properties of 
the scatter in the distance relation used. For 
example, to discuss the bias on distances es- 
timated from a TF relation given by 

I = axF^ogw + Btf (1) 

one can distinguish between the 'direct' and 
'inverse' TF relation. In equation ([^) I is the 



absolute I-band magnitude and w is the cor- 
rected HI line width. The 'direct' relation 
(called 'forward' relation by Willick 1994) is 
derived from a least-squares fit of / ClS db func- 
tion of log 1(7, while the 'inverse' relation is 
obtained from a least-squares fit of log w as a 
function of I. Independent of the distribution 
of the scatter about the line defined by equa- 
tion (|Tp, a normal distribution in logw can 
be used as an approximation for the inverse 
relation, while a normal distribution in / will 
be a better approximation for the direct re- 
lation. Although in this paper we focus our 
attention on the TF relation, our conclusions 
apply equally well to the Dn — a relation. 

A raw estimate of the distance draw from a 
measured magnitude mj and line width w 

logdraw = 0.2(m7-aTFlogw-6TF) -3 (2) 

is a biased estimate of the distance. The most 
general expression to correct draw obtained 
with the direct TF relation under the assump- 
tion that the scatter is normal in magnitudes 
is 

^ _ /o-r3n(r)exp(- ['"^-/^r^'^ )rfr 
J^r^n{r)exp{-^^^^^^^)dr 

where d^or is the corrected estimated distance, 
is the variance in the logarithm of the dis- 
tance, and n(r) is proportional to the space 
distribution of galaxies along the line of sight 
{e.g., Willick, 1994). 

Landy & Szalay (1992) have suggested a 
different expression given by 

^A n[\ndraw + ^'] , 
dcor = draw exp(— -— ), (4) 

2 h[\ndraw\ 
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where h{d) is the observed distribution of 
galaxies along the line of sight. As pointed 
out by Teerikorpi (1993), this expression is 
applicable when distances have been derived 
using the inverse TF relation. 

Both expression (|^) and can be eas- 
ily applied for the homogeneous case. They 
can also be easily applied to areas of the sky 
where the distribution of observed galaxies 
can be adequately described by one dimen- 
sional functions h{d) or n{r). However, if 
larger regions of the sky are considered, vari- 
ation with the direction in the sky must be 
taken into account. For example, estimates 
for n{r) can be obtained from redshift surveys 
taking into account expected peculiar veloci- 
ties (Hudson 1994). 

Willick (1991) also presented a formal frame- 
work to correct for selection biases which is 
similar to the corrections for Malmquist bias 
given in equation (|^). Unfortunately, his ex- 
pressions cannot directly be applied to a sam- 
ple selected as the Sc sample. 

In the analysis of the Sc sample we expect 
significant contributions from both Malmquist 
bias and selection effects due to the selection 
criteria adopted. This has prompted us to de- 
velop the MC approach presented below. In 
section 4 we compare the biases of the analyt- 
ical expressions given in this section to results 
from our MC technique discussed below. 

3. The Monte Carlo Simulation 
3.1. Method 

The ultimate goal of the method is to cre- 
ate simulated redshift - distance surveys with 
the same properties and thus the same bi- 
ases as in existing observed samples. To 
achieve this, we have to simulate properties 
hke the spatial distribution, magnitude - di- 



ameter function, color relations, distance rela- 
tions and their respective scatter properties. 
In addition, the selection criteria and com- 
pleteness of the observed samples under con- 
sideration have to be taken into account. 

For the spatial distribution of galaxies, we 
use the density field derived from redshift sur- 
veys corrected for pecuhar velocities (see be- 
low). The derived density field depends on 
the assumed f3 = Q/W'^, where Q is the den- 
sity parameter and b is the linear biasing of 
galaxies relative to matter. In this paper, we 
assume that (3 = 1. It should be emphasized 
that the density field and therefore the simu- 
lation is only slightly model dependent. 

We use this density field to create a large 
number of artificial galaxies, to which we as- 
sign properties similar to those of the ob- 
served population of galaxies. We then simu- 
late the selection of catalogs from which ob- 
servational samples are drawn. Subsequently, 
the observational sample is drawn from the 
simulated catalog, following the same selec- 
tion criteria used for existing surveys, includ- 
ing any known incompleteness of the obser- 
vational sample. Although below we assume 
that this incompleteness depends only on the 
magnitudes, this can be easily generalized to 
include incompleteness clS Qb function of direc- 
tion in the sky. Finally, estimated distances 
are derived from the assumed distance rela- 
tion. 

The final sample, which contains both the 
true and the estimated distances, can now be 
used to compute the bias field in estimated 
distance or redshift space. For each artificial 
galaxy, the bias is computed by subtracting 
from the estimated distance the known "true" 
distance of galaxies. This bias affects not only 
the position of the galaxies but also its esti- 
mated peculiar velocity which can also be cor- 
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rected using the same bias correction. Since 
the artificial catalogs of galaxies are large, the 
uncertainty in the determination of the bi- 
ases can be reduced by binning them on a 
grid either in redshift or in estimated distance 
space. This average bias can then be apphed 
to observed samples selected like the simu- 
lation. The raw estimated distance can be 
used to select the appropriate bias from the 
simulation and the estimated distance can be 
corrected for it. Peculiar velocities are sub- 
sequently computed from the bias corrected 
distances. 

A summary of the items and parameters 
that represent the built-in assumptions for 
the estimate of biases are the parameters de- 
scribing the luminosity function, the TF re- 
lation and the color-line width relation, the 
density distribution of galaxies, the peculiar 
velocities of galaxies, and the detailed selec- 
tion and completeness of the observed sam- 
ple. The main advantage of this approach 
is the fiexibility it provides to evaluate a va- 
riety of effects and their impact on the esti- 
mated distances, including the scatter in the 
various relations used, the properties of the 
distance relations and the selection and in- 
completeness of observational samples. 

3.2. Procedure 

3.2.1. Creating the 3d Galaxy Distribution 

To prepare the MC simulation, either a ho- 
mogeneous distribution of galaxies or a distri- 
bution modeled after the local universe is used 
to derive a density field in real space. For the 
latter, the Optieal-IRAS model without bias 
of Freudling et al. (1994) is used. This model 
uses the redshift distribution of optically se- 
lected galaxies in the galactic caps (6 < —30° 
or 6 > 40°), the l.QJy IRAS sample in re- 



gions |6| > 5°, and interpolates the density 
field over the galactic plane {\b\ < 5°). The 
positions of the galaxies are computed itera- 
tively from the redshifts, using the linear ap- 
proximation to estimate the peculiar velocity 
of each galaxy in each iteration. The final 
output is a self-consistent density field, which 
is stored on a Cartesian grid of 500km s^^ 
cell size. This final density field is shown by 
Willmer et al. (1995). 

The first step in the MC simulation is to 
create a list of galaxy coordinates which dis- 
tributes galaxies randomly but with density 
fluctuations identical to those in the Optical- 
IRAS model. This is accomplished by inter- 
preting the local density of the model as the 
probability p to draw a galaxy at that loca- 
tion. The X, y and z components of the posi- 
tion of potential galaxies are randomly drawn. 
The object is then rejected with a frequency 
proportional to 1 — p. 

Finally, each galaxy is given a redshift com- 
puted from its distance, the velocity field 
of the Optical-IRAS model, and an assumed 
Gaussian random component. The distances 
are assumed to be proportional to the un- 
perturbed redshift in the CMB restframe and 
are expressed in km/sec throughout this pa- 
per. For computational efficiency, the pecu- 
liar velocity vector is compiitcd only once for 
each cell on the density mesh, and all galax- 
ies within each cell are given the same pe- 
culiar velocity field. Subsequently, a Gaus- 
sian random component to the redshift of 
(Tf = 150km s~^ is added, representing an es- 
timate of the cosmic velocity dispersion. 

3.2.2. Galaxy Properties 

Most observational maps of pecuhar veloc- 
ities (e.g. Willick 1991, Courteau et al. 1993, 
or the Sc sample) are originally drawn from 
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diameter-limited catalogs such as the UGC or 
the ESO-Uppsala catalogs, and magnitude or 
other criteria are only imposed on the galax- 
ies which already have been pre-selected by 
the criteria imposed in the parent catalog. In 
order to incorporate such selection procedures 
into the MC simulation, galaxies distributed 
according to some density field must be given 
diameters, magnitudes and other parameters 
relevant for the distance relation being used. 
Sodre & Lahav (1993) have investigated the 
joint diameter-magnitude distribution func- 
tion, and have parameterized it in the follow- 
ing way: 



w.. r ^ (b - a\og(D) - B)\ 

(5) 

where B is the blue absolute magnitude 
computed from Bt, D is the diameter of the 
galaxy, and the parameters a, 6, and au 
have been derived by Sodre & Lahav. We use 
such a distribution to assign both diameters 
and blue magnitudes to the artificial galax- 
ies. For the simulation of the Sc sample, 
the parameters derived by Sodre & Lahav for 
their "S-l-I" sample with three different veloc- 
ity models were averaged. The resulting pa- 
rameters used are D* = 2231km s^^ arcmin, 
a=-4.87, b=-0.05. 

The properties of each galaxy are assigned 
as follows. First the diameters in km s~^ ar- 
cmin are drawn from an exponential distribu- 
tion function oc exp{D / D*) (see equation |^). 
For a given diameter, we derive the most 
probable value for the blue absolute magni- 
tudes B from the linear relation between di- 
ameters and magnitudes, B = a \og{D) + b, 
where a and b are the ones from the distri- 
bution functions. The final assigned magni- 
tude will later be computed by adding to B 



a Gaussian random variable with a standard 
deviation of ctm- 

Since equation (^) is given in the blue 
band, it is necessary to relate the blue Bt 
magnitude to the / magnitude. For that we 
utilize a color-line width relation of the form 
Bt — I = ac log w + be- This together with the 
I-band TF relation as given in equation (|ip 
are used to compute I and line widths from 
B. Preliminary relations have been derived 
from the Sc data and are given in table 1. As 
above, the scatter in these quantities will be 
added in the next section. 

Note that our procedure differs from that 
of Willick (1991), who assumes that the galaxy 
properties originate from a Gaussian line width 
distribution rather than a luminosity func- 
tion. Magnitudes are then computed from 
the widths via a given TF relation, and a nor- 
mally distributed scatter is added to the mag- 
nitudes. This results in a Gaussian luminosity 
function rather than the more Schechter-like 
function implicitly assumed in our approach. 

3.2.3. Correlations and Scatter 

In the previous section, the parameters Z), 
B, I and \ogw were assigned to each galaxy 
using a distribution function for D and ob- 
served relations to compute the remaining 
parameters. Up to now, we have not con- 
sidered the scatter in those relations. How- 
ever, selection biases are produced because 
galaxy samples are selected by observables 
which are closely correlated with the quan- 
tity to be measured. In particular, blue mag- 
nitudes and diameters are used later to select 
samples. Therefore, their correlation with the 
parameters of the TF relation, I band magni- 
tudes and line width, must be modeled realis- 
tically in order to obtain reliable estimates for 
the biases. The scatter in the relation used to 



6 



assign B from D equation was assumed 
to be (Ta/ = 0.709. For the TF relation, the 
scatter has been estimated to be axF = 0.3. It 
should be emphasized that the observed scat- 
ter in the TF relation and therefore the biases 
depend on the line width (e.g. Giovanelli et 
al. 1995). The value of 0.3 chosen here is the 
value measured from the Sc sample at a line 
width \ogw of about 2.4. The scatter in 
the color - line width relation, was derived 
from the Sc data to be cTc = 0.59. In table 1 
we list all the parameters used in the simula- 
tions, which are used later to determine the 
bias for samples of spiral galaxies. Finally, 
in the discussion below we also consider the 
uncertainty in the apparent blue magnitudes 
listed in the catalogs generally used to select 
the sample. Here we take this scatter to be 
am = 0.3. 

Finally, we are in the position of assigning 
values to all of the observable quantities rele- 
vant to our simulated samples such as -Bscatter, 
Scatter, and ^scatter- Thcse are computed by 
adding normal deviates to the correspond- 
ing quantities B, I and w computed before. 
For the inverse TF relation, four independent 
Gaussian random variables C,i with zero mean 
and a cr = 1 are used in the following man- 
ner. First four random variables ai to (T4 are 
defined: 

0-1 = cftf ■ ii (6) 

0"2 = O-m ■ 6 (7) 

^^3 = \joli-ol^ ■ ^3 (8) 

= \j(yl - cyl, ■ U (9) 

These random variable are added to the the 
observable quantities as computed in the pre- 
vious section in the following manner. First, 
a Gaussian with a scatter om has to be added 



to 5 so that the distribution of B and logD 
follows equation (J^). Using the random vari- 
ables defined above, we define 

-Bscatter = 5 + + CTS (10) 

in such a way that the distribution of 
-Bscatter Satisfies equation (H) and the scatter 
in the color relation is also reproduced if we 
define /scatter to be given by 

-^scatter = + + Cr4- (H) 

In order to reproduce the scatter in the TF 
relation, we define 

«^scatter=«^-10('^^+'^^+'^^)/"^^. (12) 

It can be verified that this leads to a Gaus- 
sian scatter term for B — I with a dispersion 
of (Tc by substituting ai, (72 and in the rela- 
tions. This is also approximately the scatter 
of the color - line width relation, because the 
additional terms added to w are small com- 
pared to the scatter in the color. 

Similarly, for the direct TF relation im- 
plemented by scatter in magnitudes, equa- 
tions (^) to (|T^) change to: 

^3 = V^M-^m-^TF ■ 6 (13) 

a, = ^a^^-al-^, (14) 

-Bscatter = -B + O"! + (72 + (15) 
/scatter = / + (Ti + + (T4 (16) 
W^scatter = W ■ 10('^3+.4)/a.. (^7) 
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3.2.4- Creating Simulated Surveys 

The next step in the MC simulation is to 
create galaxy catalogs similar to existing ones 
which we assume to be complete within the 
specified criteria. This step is necessary to 
simulate the exact selection procedure of the 
observer. For example, a 'magnitude-limited' 
sample drawn from the UGC catalog is not 
identical to a truly magnitude limited sam- 
ple because the inclusion of a galaxy in the 
UGC already depends on a diameter selec- 
tion threshold. Since all the observables of 
the artificial galaxies are known, this is an 
easy step. For example, for the UGC we select 
galaxies with diameters larger than 1 arcmin 
or brighter than 14.5 magnitudes. 

Next, "observed samples" of galaxies are 
drawn from the simulated catalogs. It is here 
that we can take full advantage of the MC 
approach since besides the known selection ef- 
fects (e.g. on magnitudes and sizes of galax- 
ies), even complicated selection criteria, such 
as imposed redshifts cutoffs, variations with 
direction on the sky and combination of dif- 
ferent samples can be incorporated into our 
simulations. In addition, observational sam- 
ples are characterized by an incompleteness 
function that must also be taken into account 
in building our simulated samples. In general, 
we assume that the samples are incomplete as 
a function of apparent blue magnitude only. 
The incompleteness of each sample was taken 
into account in generating our simulated sam- 
ples. 

Finally, peculiar velocities are computed 

from the estimated distances obtained from 
the TF relation and observed magnitudes, 
widths and redshifts. 



3.2.5. Computation of Bias 

The peculiar velocity of a galaxy is defined 
as Vpec — Vr — dfrue, wherc Vr iucludcs the re- 
cessional velocity and the streaming motion 
at its location, and dtme is the true velocity 
of the galaxy. This quantity is estimated from 
measurements of the redshift cz and the dis- 
1jB;I1CG draw '^pec^m — draw, where the 

measured redshift includes the cosmic random 
component crj and the estimated distance a 
large measurement uncertainty. These ran- 
dom components lead to a bias of the mea- 
sured peculiar velocity relative to Vpec- 

The final step in our MC simulation is 
to estimate these biases. This is done by 
computing raw estimates of the distances, 
corrected for homogeneous Malmquist bias, 
duMB and measured redshifts cz from the ob- 
served parameters, exactly in the same man- 
ner as in the observations. The error in the 
estimated distance 5d is computed for each 
galaxy as 

— 5d = dnUB — dfrue (18) 

Similarly, the error in the measured redshift 
relative to the underlying streaming motion 
and recessional velocity is computed as 

5z — cz — Vr- (19) 

The signs in these definitions were chosen so 
that the error in the measured peculiar veloc- 
ity is 6v ^ 5d + Sz. 

Biases are produced by the fact that the 
expectation value of these errors at a given 
estimated distance or measured redshift are 
not zero for an observed galaxy sample. In 
the MC simulation, the biases in estimated 
distances bd or in redshift space bz are com- 
puted by averaging the errors for individual 
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galaxies in the sample within spatial cubes of 
500 km s~^ size, i.e. 

hd=<Sd >cell (20) 
bz =< 5^ >cell (21) 

Depending on whether the direct or inverse 
TF relation was used, the averaging in cells 
Occii was performed differently. In case of 
the direct TF relation, galaxies were assigned 
to spatial cells based on the estimated dis- 
tances corrected for homogeneous Malmquist 
bias. This results in a bias field which con- 
tains the non-trivial biases not removed by 
the homogeneous Malmquist correction. In 
case of the inverse TF relation, this assign- 
ment to cells was made in redshift space. 

The average and root mean square of 
the biases for all N galaxies in each cell were 
estimated, homogeneous Malmquist bias cor- 
rected distances were computed and recorded. 
The resulting field of biases was then smoothed 
with a Gaussian, where the a in grid size was 
chosen to be {a^/^/N/150 km s''^)^/\ The 
minimum smoothing applied was a Gaussian 
with a cr of 1 and the maximum a used was 3. 
This results in heavy smoothing of regions 
with few galaxies and large discrepancies of 
the biases, and less smoothing in high den- 
sity regions and regions where the biases of 
the galaxies agree well. A similar smoothing 
should be applied to the real data in order to 
make the derived bias field fully apphcable. 

These bias fields are later applied to real 
observations in the following manner. First, 
the raw distance estimate of real observations 
is corrected for homogeneous Malmquist bias. 
Next, this corrected distance of the galaxy 
and its direction on the sky is used to find 
the spatial cube in which the galaxy is located 



and the corresponding biases are determined. 
The distance bias 6^ is then subtracted from 
duMB to find a final, bias corrected estimate 
of the distance. Similarly, the redshift bias hz 
is added to the measured redshift. The cor- 
rected distance and redshift are subsequently 
used to compute the peculiar velocity of the 
galaxy. 

4. Simulations 

In this section, we present a series of sim- 
ulations of "observed" peculiar velocities dis- 
tributions and compare them to the "true" 
velocities. The raw distances are computed 
from the 1-band TF relation as discussed 
above. The simulations differ in the assumed 
spatial distribution of galaxies. We first present 
a homogeneous universe and a toy model den- 
sity distribution for the purpose of demon- 
strating the general behavior of the bias and 
testing our procedure. Subsequently, we use 
a realistic density distribution, modeled after 
the density distribution derived from redshift 
surveys. Similar models will later be used 
to determine biases for existing observational 
samples. 

4.1. Tests of the MC Method 

We begin our analysis by applying our MG 
method to spatial distributions for which the 
expected biases can be computed analytically. 
This allows us to test the bias field deter- 
mined from the MG method. As our first 
model we consider a magnitude-limited sam- 
ple {niB — 14.5) drawn from a UGG-like cat- 
alog of a homogeneous distribution of galax- 
ies with no streaming motion. In this exper- 
iment, the peculiar velocities of the galaxies 
are exclusively due to the cosmic dispersion 
(7/. Figures 1 and 2 show the biases in the 
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observed peculiar velocities. In this simula- 
tion, the biases in distance 6^ are computed 
without the HMB correction to the estimated 
distances, i.e. bd =< draw - dtme >ceii- 

Figure 1 is the case corresponding to the 
use of the direct TF relation. Panel (a) shows 
the biases binned in estimated distances. Su- 
perimposed is the bias expected from a ho- 
mogeneous distribution computed from equa- 
tion (^). It can be seen that the bias behaves 
as predicted, indicating that additional selec- 
tion biases are negligible. In redshift space 
(panel b), Malmquist-type effects result from 
the cosmic random component and pro- 
duce a bias for small redshifts. The magni- 
tude of this bias is directly determined by the 
value of 150km s~^ assumed for aj. The cos- 
mic component plays the role of the error in 
distance estimate. This effect contributes to a 
negative bias because on the average the red- 
shift will be underestimated. At larger red- 
shift, this contribution becomes small and se- 
lection biases dominate. 

In figure 2 we show the bias when the in- 
verse TF relation is used. As before, the 
solid line in panel (a) represents the expected 
bias as estimated using now equation (^). 
Note that although relatively small, biases 
are present even for the inverse TF relation 
binned in redshift (panel b). At small red- 
shifts, the cause is the same as in the case 
of the direct relation discussed above. At 
large redshifts, selection biases are present 
because of the correlation between the line 
width and the blue magnitudes used to select 
the sample due to the indirect selection on line 
widths introduced by the color-line width re- 
lation. This is different from the assumption 
of Schechter (1980) that there are no selection 
effects on line widths. 

In our second model, we assume spherical 



symmetry and a radial density profile which 
simulates a void centered on the observer 
which is surrounded by a region of high den- 
sity at intermediate distances. The density 
falls off to the mean density at large distances. 
The density profile is shown in figure 3. From 
this profile, the expected streaming motion 
was computed in linear approximation using 
(3 = 1. 

In Figure 4 and 5, the measured peculiar 
velocities are plotted for the direct and for 
the inverse TF relation, respectively. In both 
figures, the measured peculiar velocities and 
their biases are plotted as a function of esti- 
mated distances and redshifts. Again, the bi- 
ases do not include the HMB correction and 
the analytical solutions computed from equa- 
tions (0) and (^ are superimposed on the MC 
results. 

Comparing figure 4 and 5 one can clearly 
see the advantage of studying the peculiar ve- 
locities in redshift space. This has previously 
been recognized by Willick (1994). Fig- 
ure 4 shows that although selection effects 
are important, they are almost independent 
of density enhancements. This is an impor- 
tant property, because in this case biases do 
not produce any apparent infall. More impor- 
tantly figure 5 shows that in the case of the 
inverse TF binned in redshift space, the biases 
are smaller than in all other combinations. 

The good agreement between the analyt- 
ical curves and both the homogeneous and 
inhomogeneous bias recovered utilizing our 
methodology demonstrates that our MC tech- 
nique is reliable. 

4.2. Realistic Model 

In order to deal with real samples, it is nec- 
essary to apply our simulations to a density 
field which adequately represents the true dis- 
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tribution of galaxies. As discussed above, we 
utilize the OpfAcal-IRAS density field for this 
purpose. In contrast to the simple cases con- 
sidered earlier, the bias bd in this section is 
computed from the HMB corrected distance 
estimates dnuB- 

The resulting velocity biases are shown in 
figure 6 as vectors along the line of sight at 
their cIhmb positions for the case of a magni- 
tude limited sample (m=14.5) and a TF re- 
lation with scatter distributed in magnitudes. 
For clarity, the scale of the bias vectors is five 
times the scale of the distances. This plot 
corresponds to a peculiar velocity field an ob- 
server would derive if there were no true pecu- 
liar velocities and the observer computed pe- 
culiar velocities from estimated distances cor- 
rected for HMB. It can be seen that the effect 
of the bias is to produce artificial infall into 
high-density regions. For instance this effect 
can be seen in the central panel near the loca- 
tions of the Great Attractor {SGX -5000, 
SGY f» -1-2000), and the Pieces- Perseus com- 
plex {SGX ^ +3000, SGY ^ -3000). 

The similarity of this map with maps of the 
predicted peculiar velocities from redshift sur- 
veys such as the one presented by Freudling et 
al. (1994) is evident. Note that these strong 
local biases average out if binned into homo- 
geneous Malmquist bias corrected estimated 
distance bins. This is shown in Figure 7. This 
shows that the HMB correction works well 
globally, but fails locally because of the large 
density variations. 

Above, we noted that the bias field in es- 
timated distance space looks similar to pecu- 
liar velocity maps predicted from the Optical- 
IRAS model reconstruction from redshift sur- 
veys. To quantify this impression, we present 
in Figure 8 the apparent velocities shown in 
figure 6 and generated by the bias versus the 



velocities predicted from the Optical-IRAS 
model. It is such a plot which in principle can 
be used to derive a value for f3. The compari- 
son of the distribution of points with the lines 
for different /3 suggests that from the data one 
would derive a /3 in the interval 0.2-0.5 even 
when the true peculiar velocities are zero, i.e. 
the true (3^0. This demonstrates that f3 
cannot be derived from data which have been 
corrected only for HMB. 

In the previous section, it was argued that 
distances estimated with the inverse TF re- 
lation are not dominated by Malmquist bias 
if they are plotted function of redshift. 
Instead, such distances are biased by selec- 
tion effects (see Figure 2). In order to demon- 
strate this difference, we plot in Figure 9 the 
biases for a model identical to the one shown 
in Figure 6 but with the scatter in the TF re- 
lation taken to be in the line width, and the 
biases are shown in redshift space. It can be 
seen that the biases show a completely differ- 
ent behavior than those in estimated distance 
space shown in Figure 6. They are mainly a 
function of distance, and do not concentrate 
around high density regions. For the chosen 
magnitude limited sample of 14.5, the biases 
are significantly smaller than the one seen in 
figure 6. 

5. Applications 

While the purpose of the previous dis- 
cussion was to understand the properties of 
Malmquist bias in combination with selection 
effects, the purpose of this section is to pre- 
dict the magnitude of these effects for real sur- 
veys. Therefore we create simulated redshift- 
distance surveys which correspond to samples 
actually observed. A summary of the sur- 
veys considered and their selection criteria is 
listed in table 2. The samples considered are 
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those of Lynden-Bell et al. (1988, hereafter 
7S), Wilhck (1990, hereafter W90) and the 
new Sc sample. 

5.1. 7S Sample 

The parameters for the simulation of the 
7S sample of elliptical galaxies are given in 
table 3 and the incompleteness func- 
tion of magnitude was taken from Faber et al. 
(1989). Since the TF relation is substituted 
by a Dn — a relation (see table 2), the color- 
line width relation must also be replaced by 
the equivalent relation. We could use a re- 
lation between D — Dn and o", derived from 
the data, which plays the same role as the 
color-line width relation and allows for possi- 
ble correlation between the parameters used 
to select the sample with both of the quan- 
tities of the distance relation. No significant 
correlation between D — Dn and a was found, 
therefore Dn was computed from D by adding 
the constant and Gaussian scatter given in ta- 
ble 3. 

In figures 10, the derived bias field for 
the 7S sample is shown. The biases are 
again shown in estimated HMB corrected dis- 
tances space. The biases produce some ap- 
parent infall into the 'Great Attr actor' (at 
SGX ^ -4200km s-\ SGY ^ +800km s-^ 
SGZ ^ —700km s~^), but the effect is not 
much stronger than in other parts of the sky. 
In order to address the question whether such 
biases could produce the velocity field as sug- 
gested by 7S, we bin in figure 11 the biases 
as a function of the cosine of the angular dis- 
tance from the GA in several distance shells. 
The binning and presentation corresponds ex- 
actly to the one used by 7S in their figure 5. 
It can be seen that the biases binned in such 
a manner are independent of the angular dis- 
tance from the GA, unlike the behavior the 7S 



find in their data. We conclude that residual 
biases do not significantly affect the conclu- 
sion of 7S regarding the existence of a large 
infall towards the GA within the volume con- 
sidered. 

5.2. WOO Sample 

In order to estimate the incompleteness of 
the W90 sample, we have binned the data in 
0.5 magnitudes intervals and compared the re- 
sulting counts to those obtained in the CGCG 
(Zwicky et al. 1961-68) catalog. This incom- 
pleteness was imposed on the simulated sam- 
ple. Subsequently, the biases were computed 
using the direct TF relation as discussed in 
section 3.3. It should be noted that the TF 
relation we use differs from that assumed by 
Willick (1990). In principle we could use the 
same quadratic relation, but we believe that 
the major source of uncertainty stems instead 
from the lack of detailed information regard- 
ing the properties of the scatter in the TF 
relation. In particular, the sample contains 
galaxies as faint as magnitude 16.5, but the 
number of galaxies in each magnitude inter- 
val is small. This leaves the possibility of sub- 
stantial unrecognized systematic effects in the 
incompleteness of the sample. 

In figure 12, the biases for the simulated 
W90 sample is shown. The figure shows the 
biases both function of redshift and as 
a function of (Ihmb- The biases are shown 
separately for models with no true peculiar 
velocities and with a velocity field given by 
the Optical-IRAS model. The biases are es- 
sentially independent of the assumed velocity 
model. A comparison of the two panels in fig- 
ure 12 shows that for a sample selected hke 
W90, the biases out to 5000km are smaller 
if the peculiar velocities are plotted func- 
tion of duMB rather than redshifts. Beyond 
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that distance, the bias is similar in amphtude 
but with opposite sign. 

In figure 13, the W90 data are corrected 
for the bias shown in figure 12 and binned in 
bins of 500 km s^^. Panel (a) corresponds 
to figure 1 shown by W90, whereas panels (b) 
and (c) show the bias corrected data as a func- 
tion of redshift and distance, respectively. A 
significant infall into the Pisces-Perseus su- 
percluster complex can be seen both in pan- 
els (b) and (c). While the corrections (as 
shown in figure 12) are substantial, the basic 
conclusion of W90 that there is a large bulk 
motion of the whole region in the direction 
of the Local Group cannot be ascribed to the 
biases discussed in this paper. 

5.3. Sc Sample 

Samples selected using a redshift criterion 
such as the Sc sample are a completely differ- 
ent case. The Sc sample also differs from the 
previously discussed sample by the fact that 
a large fraction of all galaxies which satisfy 
the selection criteria have actually been ob- 
served. This lowers the probability for large 
biases due to unrecognized systematic effects 
in the selection of the sample. In figure 14 the 
'measured' peculiar velocities of a sample of 
galaxies selected with the Sc sample criteria is 
shown function of estimated distance and 
redshift. As before, the estimated distances 
in this this figure have been HMB corrected. 
At small distances (up to about 6000km/sec) 
the biases in addition to the HMB are small. 
This refiects the changing diameter limit as 
a function of redshift adopted in selecting the 
sample. On the other hand, at large distance 
the peculiar velocities are strongly affected by 
the redshift cutoff which leads to the depopu- 
lation of galaxies with positive peculiar veloc- 
ities. The full bias map which we will actually 



use in our analysis of the Sc sample is shown 
in figure 15, again plotted against HMB cor- 
rected distances. 

6. Summary 

We have presented a new approach to cor- 
rect observed peculiar velocity fields. The 
Monte-Carlo approach we have adopted al- 
lows us to simultaneously include all known 
sources of biases in a realistic way. The result- 
ing simulated biases can be applied to exist- 
ing observed samples with complex but well- 
defined selection criteria, such as the new Sc 
sample. 

We have tested the MC simulation in sim- 
ple cases, for which bias correction schemes 
are known, and recovered the expected biases 
with high accuracy. We have also applied our 
technique to real samples and computed the 
expected biases. We found that the results 
presented by Lynden-Bell et al. are virtually 
free of biases. On the other hand, the data as 
presented by Willick (1990) are significantly 
biased, but the biases do not affect his main 
conclusion. Finally, we derived the expected 
bias for the new Sc sample of Haynes et al. 
(1993) that will be used in future analysis of 
this sample. 
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Figure Captions 

Figure 1: Biases of an homogeneous distri- 
bution of galaxies observed witli 14.5 magni- 
tude limit drawn from the UGC catalog. The 
scatter of the TP relation of this sample is 
distributed in m, corresponding to the direct 
TP relation. The velocities were binned in 
500km s~^ bins in estimated distance d^^w in 
panel (a) and redshift cz in panel (b). The 
units for distances, redshifts and biases are 
km/sec. The solid line in the panel (a) is the 
homogeneous Malmquist bias computed from 
equation (||). 

Figure 2: Same figure as Figure 1, but with 
the scatter in the TF relation distributed in 
logw, corresponding to an inverse TF rela- 
tion. The solid curve in panel (a) is the ex- 
pected bias for the observed distribution of 
galaxies computed from equation (|). 

Figure 3: Density profile and peculiar veloc- 
ities of the toy model described in the text. 
Figure 4: Observed peculiar velocities (pan- 
els a and c) and biases (panels b and d) 
from the toy model for the direct TF rela- 
tion. Panels (a) and (b) show the velocities 
and biases binned in estimated distances draw, 
whereas panels (c) and (d) show the them 
binned in redshift. The dashed lines in pan- 
els (a) and (b) are the expected homogeneous 
Malmquist bias. The solid curve in panel (b) 
is the bias as computed from equation (|^), 
using the known density profile as shown in 
figure 3. Note the strong bias for the redshift 
binning (panel d). 

Figure 5: Observed peculiar velocities and 
biases from the toy model for the inverse TP 
relation. The presentation is the same is in 
figure 4. The solid curve in panel (b) is the 
bias as computed from equation (^), using the 
observed distribution of galaxies derived from 



the data. Note that there is some bias even 
for the redshift binning (panel d). 

Figure 6: Biases of a 14.5 magnitude lim- 
ited sample shown function of HMB cor- 
rected positions. The biases were computed 
using the direct TP relation and are shown 
as vectors at their positions in Supergalactic 
Cartesian coordinates. This plot corresponds 
to a peculiar velocity field an observer would 
derive if there were no true peculiar velocities 
and the observer computed peculiar velocities 
from estimated distances corrected for HMB. 
The units for distances are km/sec. The scale 
of the velocity vectors is five times the scale 
of distances. 

Figure 7: The biases shown in figure 6 
binned in HMB corrected estimated distances 
bins. Note the very small amplitude of the 
bias if averaged over the whole sky. 

Figure 8: Estimated HMB corrected pecu- 
liar velocities for a model without true pecu- 
liar velocities plotted versus the peculiar ve- 
locity as predicted by the Optical- IRAS model 
for (3 = 1. If no biases were present, the two 
quantities should not correlate. The observed 
correlation between estimated and predicted 
velocities could be mistaken as evidence for 
gravitationally induced streaming motion if 
biases are not taken into account. The lines 
show expected correlations for several differ- 
ent values of j3. The observed correlation is 
similar to the one expected for values of (3 
between 0.2 and 0.5. 

Figure 9: Biases of distances derived from 
the TP relation for a 14.5 magnitude limited 
sample shown in redshift space. The bias map 
shown here is again plotted function of 
Supergalactic coordinates. The velocity vec- 
tors are five times the scale of distances. This 
is to be compared to figure 6, which shows a 
similar plot shown in corrected estimated dis- 
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tance space, which has a completely different 
behavior of the biases. 

Figure 10: Simulated biases for the 7S sam- 
ple, binned in estimated distance space. The 
scale of the arrows is 5 times the scale of dis- 
tances. 

Figure 11: The unsmoothed biases from fig- 
ure 10 plotted as a function of the cosine of 
the angular distance from the Great Attrac- 
tor in 3 estimated distance ranges, similar to 
figure 5 in 7S. It can be seen that conclusions 
by 7S based on this plot are not subject to 
biases. 

Figure 12: Biases computed for the simu- 
lated sample of Willick (1990). The panel (a) 
shows the biases as a function of redshift. 
Panel (b) shows the biases as a function of 
HMB corrected distances. The error bars 
show the biases computed for a model with 
a velocity field as described in the text. The 
solid points show the biases for a model with- 
out any 'true' velocities. The sohd points are 
offset horizontally from their true position by 
100 km s~^ in order to make them easier to 
distinguish from the error bars. 

Figure 13: The W90 sample corrected for 
biases and binned in bins of 500 km s~^. 
Panel (a) shows the raw data as a function 
of redshift in the CMB restframe. Panel (b) 
shows the corrected peculiar velocities binned 
in redshift bins of 500km s^^ size. Panel (c) 
shows the corrected data binned in equivalent 
estimated distance bins. 

Figure 14: HMB corrected estimates of the 
peculiar velocities of a sample of galaxies se- 
lected like the Sc sample, binned in estimated 
distances (panel a) and redshifts (panel b). 
The small points indicate individual galax- 
ies, whereas the large points are averages in 
bins of 500km s~^. The dashed fine is the 
computed homogeneous Malmquist bias. The 



most prominent feature when plotted as a 
function of the estimated distance is the red- 
shift cutoff at 7500 km s^-*^, which has the 
effect that galaxies with a redshift > 7500 
km s~^ are only included if they have neg- 
ative peculiar velocities. 
Figure 15: Bias map for the simulated sam- 
ple of Haynes et al. (1994). The unit of the 
velocity vectors is identical to the unit of dis- 
tances. The vectors are plotted at the posi- 
tion of the estimated distances corrected for 
homogeneous Malmquist bias. 
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Table 2: Selection Criteria of Samples 



Sample catalogs type area mag. sizes redshifts 

~~Sc UGC Sc \b\ > 10° 5' >d> 2.5' < 3000 

ESO 5' > d > 1.5' 3000 to 5000 

MGC 5' > d > 1.3' 5000 to 7500 

7S RC2 E \b\ > 2.5° < 13.0 d>2' < 8000 

RSA 

UGC 

ESO 

Willick UGC Spirals 70 < Z < 160° < 16.5 d > 0.6' 
CGCG b < 
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Table 1: Model Parameters for Spiral galaxies 



Luminosity 


function 


D* km s"^ 


2231 


a 


-4.87 


b 


-0.05 




0.709 


TF 


relation 


(Ttf 


0.3 


O'TF 


-7.48 


brF 


-2.53 




0.3 




1.78 


be 


-2.54 




0.59 



Tabic 3: Model Parameters for E + SO 



Luminosity 


function 


D* km s"^ 


2489 


a 


-4.62 


b 


-1.66 




0.396 


Dn-cr 


relation 




0.6 


OjDn—a 


1.2 




0.3 


D-D^ 


0.05 




0.5 



Figure 1: Biases of an liomogeneous distribution of galaxies observed with 14.5 magnitude limit 
drawn from the UGC catalog. The scatter of the TF relation of this sample is distributed in m, 
corresponding to the direct TF relation. The velocities were binned in 500km s~^ bins in estimated 
distance draw in panel (a) and redshift cz in panel (b). The units for distances, redshifts and biases 
are km/sec. The solid line in the panel (a) is the homogeneous Malmquist bias computed from 
equation (^. 



Figure 2: Same figure as Figure 1, but with the scatter in the TF relation distributed in \ogw, 
corresponding to an inverse TF relation. The solid curve in panel (a) is the expected bias for the 
observed distribution of galaxies computed from equation (H). 



Figure 3: Density profile and pecufiar velocities of the toy model described in the text. 



Figure 4: Observed peculiar velocities (panels a and c) and biases (panels b and d) from the 
toy model for the direct TF relation. Panels (a) and (b) show the velocities and biases binned 
in estimated distances draw) whereas panels (c) and (d) show the them binned in redshift. The 
dashed lines in panels (a) and (b) are the expected homogeneous Malmquist bias. The solid curve 
in panel (b) is the bias as computed from equation (^, using the known density profile as shown 
in figure 3. Note the strong bias for the redshift binning (panel d). 



Figure 5: Observed peculiar velocities and biases from the toy model for the inverse TF relation. 
The presentation is the same is in figure 4. The solid curve in panel (b) is the bias as computed 
from equation (H), using the observed distribution of galaxies derived from the data. Note that 
there is some bias even for the redshift binning (panel d). 



Figure 6: Biases of a 14.5 magnitude limited sample shown as a function of HMB corrected 
positions. The biases were computed using the direct TF relation and are shown as vectors at 
their positions in Supergalactic Cartesian coordinates. This plot corresponds to a peculiar velocity 
field an observer would derive if there were no true peculiar velocities and the observer computed 
pecuhar velocities from estimated distances corrected for HMB. The units for distances are km/sec. 
The scale of the velocity vectors is five times the scale of distances. 



Figure 7: The biases shown in figure 6 binned in HMB corrected estimated distances bins. Note 
the very small amplitude of the bias if averaged over the whole sky. 



Figure 8: Estimated HMB corrected peculiar velocities for a model without true peculiar velocities 
plotted versus the peculiar velocity as predicted by the Optical-IRAS model for /3 = 1. If no biases 
were present, the two quantities should not correlate. The observed correlation between estimated 
and predicted velocities could be mistaken as evidence for gravitationally induced streaming motion 
if biases are not taken into account. The lines show expected correlations for several different values 
of p. The observed correlation is similar to the one expected for values of /3 between 0.2 and 0.5. 



Figure 9: Biases of distances derived from the TF relation for a 14.5 magnitude limited sample 
shown in redshift space. The bias map shown here is again plotted as a function of Supergalactic 
coordinates. The velocity vectors are five times the scale of distances. This is to be compared 
to figure 6, which shows a similar plot shown in corrected estimated distance space, which has a 
completely different behavior of the biases. 



Figure 10: Simulated biases for the 7S sample, binned in estimated distance space. The scale of 
the arrows is 5 times the scale of distances. 



Figure 11: The unsmoothed biases from figure 10 plotted as a function of the cosine of the 
angular distance from the Great Attractor in 3 estimated distance ranges, similar to figure 5 in 
7S. It can be seen that conclusions by 7S based on this plot are not subject to biases. 



Figure 12: Biases computed for the simulated sample of Willick (1990). The panel (a) shows 
the biases as a function of redshift. Panel (b) shows the biases as a function of HMB corrected 
distances. The error bars show the biases computed for a model with a velocity field as described 
in the text. The solid points show the biases for a model without any 'true' velocities. The solid 
points are offset horizontally from their true position by 100 km s~^ in order to make them easier 
to distinguish from the error bars. 



Figure 13: The W90 sample corrected for biases and binned in bins of 500 km s^^. Panel (a) 
shows the raw data as a function of redshift in the CMB restframe. Panel (b) shows the corrected 
pecuhar velocities binned in redshift bins of 500km s"^ size. Panel (c) shows the corrected data 
binned in equivalent estimated distance bins. 



Figure 14: HMB corrected estimates of the peculiar velocities of a sample of galaxies selected 
like the Sc sample, binned in estimated distances (panel a) and redshifts (panel b). The small 
points indicate individual galaxies, whereas the large points are averages in bins of 500km s~^. 
The dashed line is the computed homogeneous Malmquist bias. The most prominent feature when 
plotted as a function of the estimated distance is the redshift cutoff at 7500 km s~^, which has the 
effect that galaxies with a redshift > 7500 km s~^ are only included if they have negative pecuhar 
velocities. 



Figure 15: Bias map for the simulated sample of Haynes et al. (1994). The unit of the velocity 
vectors is identical to the unit of distances. The vectors are plotted at the position of the estimated 
distances corrected for homogeneous Malmquist bias. 



